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ABSTRACT 



Various early-type dwarf galaxies with disk features have been identified in the Virgo cluster, including objects that display weak 
grand-design spiral arms despite being devoid of gas. Are these still related to the classical dEs, or are they a continuation of ordinary 
spiral galaxies? Kinematical information of acceptable quality is available for one of these galaxies, IC 3328. We investigated its 
dynamical configuration, taking the effect of asymmetric drift into account, and using the Toomre parameter, as well as density wave 
considerations. The derived mass-to-light ratios and rotational velocities indicate a significant dynamically hot component in addition 
to the disk. However, any unambiguous conclusions will need to await further data for this and other early-type dwarfs with spiral 
structure. 



Key words. Galaxies: dwarf - Galaxies: elliptical and lenticular, cD - Galaxies: fundamental parameters - 
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1. Introduction 



Galaxies: kinematics and 



As the most numerous type of galaxy in clusters and the pos- 
sible descendants of the building blocks in hierarchical struc- 
ture formation, early-type dwarf (dE) galaxies play a key role 
toward understanding galaxy cluster evolution. Initially believed 
to be spheroidal objects having old stellar populations and pre- 
ferring the high-density regions of clusters, today's picture of 
dEs is far more diverse. In the Virgo cluster, several subclasses 
with significantly different characteristics exist (Lisker et al. 
2007), and these are correlated with environmental density. 
Those dEs that populate less dense cluster regions par tly have 
younge r stellar populations (Lisker et al. 2006a, 2008), flatte r 
shapes dFereuson & Sandagdll989T lLisker et all l2006bl l2007l) . 
and clustering properties not like giant ellipticals but similar to 
spiral galaxies dLisker et"aT1l2007l) . 

Among them is the "dE(di)" subclass, characterised by weak 
disk features - like spiral arms or bars - that could only be seen 
through unsh arp masks or by sub tracting a model of the smooth 
galaxy light (Lisk er et al.l l2006bb . The first discovery of spiral 
structure in a dE was reported by iJerien et all d2000h for the 
galaxy IC 3328 (VCC 856), which is also the focus of the work 
presented here (Fig.Q]). The dE(di)s are not a negligible species, 
but they make up one third of the brighter (Mb < -16) Virgo 
cluster dEs, reaching a fraction of 50% at the bright end. For 
those showing spiral substructure, the arm opening angles are 
inconsistent wit h being the mere remainders of late -type pro- 
genitor galaxies ( Jerien et al. 2000l ILTsker et al.ll2006bl) . and the 
arms do not show up i n colour maps that w ould indicate stellar 
population differences (Lisker et al. 20 06bl) . 




* Partly based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile, for pro- 
gram 077.B-0785. 
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Fig. 1. Appearance and light profile of IC3328. Surface 
brightness in SDSS-r versus radius, a dopting a distance mo dulus 
m-M = 3 1.0 mag (d = 15.85 Mpc. Ivan den Berghll 19961) . cor- 
responding to a scale of 77 pc/". The half-light radius r e , as well 
as 2.2 scalelengths h (using r e = 1.68/z), are indicated by the 
vertical dotted lines. The images shown on the right are, from 
top to bottom, the SDSS r-band image, a deep white-filter im- 
age obtained with ESO 2.2m/WFI (program 077.B-0785), and 
an unsharp mask image of the latter, revealing the spiral arm 
structure. The half-light aperture is indicated in the SDSS and 
unsharp mask image. 



An indication that the dE(di)s might have been formed out 
of infalling disk galaxies is that spiral arms and bar struc- 
tures in d Es also arise in N-bod y simulations of galaxy ha- 
rassment (Mastropi etro et ail 120051) . in which late-type galax- 
ies accreted by the cluster experience a violent structural trans- 
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Fig. 2. Light profile in z-direction of IC3435. SDSS r-band 
image (inset) of IC 3435, and the intensity profile obtained by 
collapsing the outlined box-shaped region along the galaxy's 
major axis. The box has a side length corresponding to the ma- 
jor axis of the half-light elliptical isophote, or equivalently, three 
times its minor axis. Smearing due to seeing effects can occur in 
the innermost ~ ±1.5". 



formation dMoore et al] 1 19961: also see iKormendy et al.1 12009). 
Alternatively, the dE(di)s - being compact, gas-poor disk galax- 
ies - might simply constitute the low-luminosity counterpart to 
normal SO/Sa galaxies, with a few "dwarf-like" SO/Sa galaxies 
possibly bridging the gap to the more luminous systems (Lisker 
etal. 2006a). 

Do the dE(di)s have a significant dynamically hot compo- 
nent, which would be expected for dEs with embedded disks, 
or are they consistent with being pure disk galaxies? A first 
assessment can be provided by examin ing an apparently e dge- 
on dE(di), IC3435 (VCC 1304; Fig. l2t ILiske7etani2006bh . Its 
axis ratio of 0.33 e quals the estimate fo r the intrinsic axis ratio 
of the dE(di)s from Lis ker et alj d2007l) . While it is, of course, 
not known whether this galaxy hosts similar spiral structure to 
IC 3328, it still appears worth examining whether a thin disk 
component is visible in the vertical light profile, i.e. perpendicu- 
lar to the disk. However, from the SDSS r-band image, no such 
component is seen (Fig.[2j, and the profile appears perfectly con- 
sistent with a single exponential decline. 

In this context, it would be obvious to ask about the kine- 
matical properties of the dE(di)s with spiral arms. Unfortunately, 
useful kinematical data is only available for one of these enig- 
matic objects, namely IC 3328, and only along its major axis 
(Simi en & Prugniell2002l) . Here we attempt to draw conclusions 
on the nature of IC 3328 from these data. 



2. Data 

2. 1 . Structural measurements 

Total magnitude, half-light radius, and ellipticity of IC 3328 
were measured from r-band images of the Sloan Digital 



Sky Survey dAdelman-M cCarth v et alj I2007I) . as outlined in 
ILisker et all d2008l) . The ellipticity of 0.87, measured at the half- 
light radius, is consistent with the ellipticity that could be derived 
from the spiral structure only, as illustrated in Fig.Q] This yields 
an inclination i = 30°. The radial surface brightness profile in 
r, shown in Fig.Q] was meas ured from annuli of fixed elliptical 
shape, using IRAF/ellipsf] (Todv 1993). Out to the radial extent 
of the rotation curve (Fig. [3]) it would be compatible with an ex- 
ponential profile, i.e. a straight line. Nevertheless, a change in 
the profile slope is seen at larger radii. This is included in the 
discussion (Sect. [4]). 

2.2. Stellar mass-to-light ratio 

From the analysis of S PSS multicolour photometry for Virgo 
dEs (ILisker et al.ll2008l) . we can derive an estimate of the stellar 
mass-to-light ratio (M/L) of IC3328, using the g — r and i - z 
colours measured within th e half-light radius. With st ellar pop- 
ulation synthesis models of iBruzual & Chariot! d2003l). ' 'Padova 
1994" isochrones dBertelli et al.lll994l) . a Chabrierl d2003h initial 
mass function, and an exponentially declining star formation rate 
(decay time r = 1 Gyr), we derive luminosity-weighted age and 
metallicity values of 7.5 Gyr (time since the onset of star forma- 
tion) and Z = 0.004 (corresponding to [Fe/H] = -0.64). When 
using simple stellar population (SSP) models instead, the age 
changes to 5.9 Gyr. This perfectly agre es with the spectroscopi c 
results from the Lick index analysis of Michielsen et al. (2008), 
who find 6.0 Gyr and [Fe/H] = -0.64 using the SSP models of 
Vazd ekis et al.1 d!996l) . It also agrees with a preliminary spectro- 
scopic analysis of ESO VLT/FORS2 spectra (Paudel et al., in 
prep.). 

The corresponding stellar M/L in V is (M/L)y = 1.4. 
Despite the good agreement of age and metallicity from differ- 
ent sources, we use the rather large age uncertainties given by 
iMichielsen et alj d2008l) (6+4 Gyr) to obtain M/L uncertainties 
(1.4+Qg). This conservative approach is also supposed to cover 
any potential M/L uncertainties inherent to the modelling pro- 
cess. 



2.3. Kinematical data 

A rotatio n and velocity disp e rsion curve of IC 3328 was pub- 
lished by Simien & Prugniel (2002), obtained from major-axis 
long-slit spectroscopy taken at the Observatoire de Haute- 
Provence. For the analysis presented here, the average of pub- 
lished velocity data points at opposite sides of the centre was 
taken, and errors propagated accordingly. We correct for the in- 
clination (Sect. 12. It and consider the effect of finite disk thick- 
ness on the line-of-sight velocity distribution (Appendix lAl . 
yielding the rotation curve shown in Fig. [3] With the adopted 
intrinsic axis ratio of 0.33 dLisker et al.ll2007l) . the resulting ve- 
locity correction is negligible (< 3%) for ra dii beyond . 5 disk 
scalelengths. The rotation curve published bv lGeha et al.1 (2003) 
is not used due to its significantly smaller radial extent. 



3. Kinematical analysis 

Our approach is to analyse the kinematical data of IC 3328 with 
the assumption of a pure disk, as implied by the vertical light 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 3. Rotation curve and the effect of asymmetric drift. Large filled circles are the inclination-corrected measurements from 
Simien & Prugniel (2002). The lower solid curve is a fit of f(r) = a(l - e~J) to them. Taking the asymmetric drift into account (see 
text) results in the small filled circles, along with the corresponding fit using the same function (upper solid c urve). In the left pa nel, 
we use <x R = 45 km/s, while 71 km/s is used in the right panel. The inclination-corrected measurements from lGeha et al.l d2003l) are 
shown as filled squares in the left panel. 



profile of IC 3435 (Sect. [TJ. We then discuss whether a two- 
component galaxy would be more likely based on the results. 

3.1. Velocity dispersion 

The observed line-of-sight velocity dispersion cn . s . is more or 
less constant with radius, within the measurement uncertainties 
( Simi en & Prug niel 2002), at 35 km/s. For major axis spectra, it 
can be expressed in terms of the angular and the vertical velocity 
dispersion as 

(1) 



'Las. 



cr fI . sin i + o\ cos ; 



Following iGerssenl (l2000h . we adopt cr z = 0.8<x R as our "work- 
ing value", but also perform our analysis using factors of 0.4 
and 0.9 instead. This range covers the values found for galaxies 
of type T=l to 5. For cr , we have 



1 + - dVc 

v c dR 



(2) 



CTj, _ 1 

OR V2 

which we estimate to be between 0.7 and 1 . With the above range 
of factors, this yields values for cn s between 0.5 and 0.9 times 
<x R , leading to 38 < cr R < 71 km/s, with our working value being 
45 km/s. 



3.2. Asymmetric drift 

While the asymmetric drift has a minor effect on the observed 
rotation curves of normal spiral galaxies, it can be much more 
significant for dwarf galaxies. The true orbital velocity at a given 
radius, v c (R), and the observed velocity, i>$, are related b y the 
radial Jeans equation as follows (Binnev & Tremaine 1987): 

R-— — ^— + 1 + — % ^ = (3) 

dR at 



with v c - the asymmetric drift velocity. (We neglect here the 
tilt of the velocity ellipsoid.) Adopting da R /dR = 0, this leads 
to 



1 / R dv c \ Rdl.\ 
2\ ~ V c ~dRj + XdR) 



(4) 



with the surface mass density E. 

Since we do not want to impose any kinematical model 
on the analysis, we adopt the following simple method to de- 
rive dv c /dR and determine v c (R). The observed, inclination- 
corrected rotation curve is fitted with f(r) = a(l - e~*), and 
its derivative is used as initial estimate for dv c /dR at each data 
point. The thus derived first iteration for v c (R) is then also fitted 
with f(r), and the derivative is now used as a better estimate for 
dv c /dR, applying equation [4] Here, the surface mass density E 
is taken to be the observed luminosity density, which is obtained 
from the surface brightness profile (Fig.[T|) and the adopted dis- 
tance of 15.85 Mpc. No specific M/L value needs to be assumed 
for the conversion, since it occurs in both the enumerator and the 
denominator of the last term of equation |H thus cancelling out. 
In this and all following steps, we exclude data points at R < 3", 
since the luminosity density is only reliable for radii significantly 
larger than the SDSS seeing FWHM (» 1 .5"). 

This leads to the second iteration on v c (R), and is repeated 
one more time. The resulting data points for v c (R) are then also 
fitted with /(r), yielding a rotation curve instead of discrete data 
points alone. The differences between the last and second-last 
iterations are very small (< 2 km/s). In this procedure we always 
set u c (0) = 0. The resulting curves are shown in the left panel 
of Fig. [3]for our working value <r R = 45 km/s (Sect. 13.11 ), and in 
the right panel for the extreme case of <x R = 71 km/s. These are 
the rotation curves that enter the calculations of the following 
sections. 
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Table 1. Mass-to-light ratio i n V, u sing the colour transforma- 
tion V = r + 0.27 (Smi th et alj|2002|) . 



Case 


M/L v 

(<r R =45km/s) 


M/Ly 

(a* =71 km/s) 


(°T{,Alt=10km/s) 


2=1 


11.2 


25.0 


0.9 


2 = 1.5 


7.4 


16.7 


0.6 


2 = 2 


5.6 


12.5 


0.4 


« = ^max 


2.9 


6.2 


0.8 


Stellar 




1 4 +L4 





The alternative scenario, in which only part of the observed velocity 
dispersion is attributed to the disk (see Sect.[4]l. 



Table 2. Rotational velocity, with values at the half-light radius 
for an exponential disk and a Mestel disk, using the same cases 
as in Table Q] 



Case 


u c 1 km s 1 


v c /km s 1 


w c /km s 1 a 




(tr R =45 km/s) 


(ctr =71 km/8) 


(°"R,AIt=l 0km / s ) 


Exp., Q = 1 


125 


188 


38 


Exp., 2 = 1-5 


102 


153 


31 


Exp., Q = 2 


89 


133 


27 


Exp., A = A mllx 


64 


94 


36 


Measured + asym. drift 


59 


87 


28 


Mestel, Q = 1 


113 


169 


35 


Mestel, 2 = 1-5 


92 


138 


29 


Mestel, g = 2 


80 


119 


25 


Mestel, A = A max 


58 


84 


34 



The alternative scenario, in which only part of the observed velocity 
dispersion is attributed to the disk (see Sect.f4]l. 



3.3. Toomre stability parameter 

An important d iagnostic of the dynamics of a galactic disk is the 
iToomrd 0964) stability parameter: 



Q = 



3.36GI 



(5) 



where S is the surface density and k the epicyclic frequency. The 
value of Q for a thin disk must lie in the range 1 < Q < 2, 
to ensure that the disk is not only stable against l ocal Jeans in - 
stabilities, but also able to develop spiral arms (lFuchsll2Q01h . 
By choosing a value for Q, the mass density E can be obtained, 
yielding a mass-to-light ratio when combined with the observed 
luminosity density. 

To properly consider the disk thickness (adopted intrinsic 
axis ratio 0.33, see Sects. Q]and[2]), we calculate how much more 
mass a thicker disk contains per surface area if it has the de- 
sired properties for spiral arm development (Appendix IB1. This 
is equivalent to using an "effective Q" value that is lower than 
the actual one for a thin disk. The above range of 1 < Q < 2 thus 
changes to 0.62 < Q e g < 1 .24 for the case cr R = 45 km/s, and 
0.64 < Q eff < 1.28 for cr R = 71 km/s. 

Table [T] shows the resulting values for M/L at the half- 
light radius when choosing Q = 1.0,1.5,2.0 (i.e., <2eff = 
0.62,0.93, 1.24 for <x R = 45 km/s, etc.). These are compared to 
the stellar M/L from Sect.|2~!2l 



The thus derived mass density at the half-light radius can be 
used to predict the rotational velocity for an exponential disk 
model and for a Mestel disk. These values are listed in Table [2] 
and illustrated in Fig. [4] again for the three Q values. 



3.4. Density wave theory 

An alternative estimate for the mass density can be obtained 
from density wave theory. It predicts the circumferential wave- 
length A max at whic h density wave s mainly grow. In terms of the 
critical wavelength dToomrdl 19641) . 



^nrit — 



4tt 2 G2 



(6) 



the preferred wavelen gth is given by dToomrd Il981h 
lAthanassoula et al.| [T987) 



'Imax — X I 77" 
I "0 



(7) 



The coefficient X depends on the shape of the rotation curve, 
measured by Oort's constant A, with 



A _ 1 / R dv c 
~ 2 1 ~ V c ~dR 



(8) 



and is given in iFuchsl (HbOl)EI The number of spiral arms, m, is 
given by how often /i max fits onto the annulus (Fuchs 2008): 



2nR „ 
m = = 2 



From equations [6] [7] and [9] we have 
1= K2R 



4 - G -*(i) 



(9) 



(10) 



In contrast to the Toomre parameter, this involves no veloc- 
ity dispersionQ An approximation for taking the actual thickness 
of the disk into account yields a 15% lower value for the coeffi- 
cient X (see AppendixO. which we use in the calculation of the 
resulting M/L values at the half-light radius (Table Q] labelled 
"A = /Imax")- Like above for the Toomre parameter, rotational 
velocities can be predicted for an exponential and a Mestel disk 
(Tableland Fig. [4]!. 

4. Discussion 



The analysis of IC 3328 by Uerien et all d2000t) shows that 
the spiral structure holds clues to its dynamical configuration. 
However, no rotation curve had be en obtained then. With th e 
now available kinematical data from Simien & Prugniel (2002), 
a more detailed investigation was possible. 

While the mass-to-light ratios derived with the Toomre cri- 
terion clearly exceed the stellar M/L (Table [TJ, the value de- 
rived using density wave theory is still compatible with the stel- 
lar value. The predicted rotation curve of the exponential model 
from A = /l max provides good agreement with the measured, 
asymmetric-drift-corrected curve (Fig. |4|. However, as for the 



2 Actually, l/X is given there. 

3 Since these calculations are performed based on the asymmetric- 
drift-corrected rotation curve, a dependence on velocity dispersion is 
implicitly included through k, since the derivation of l> c from the ob- 
served Do involves <tr (see equation!?]). 
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Fig. 4. Measured versus predicted rotational velocities. The 

black lines are the same as in the left panel of Fig. [3] i.e. for 
cr R = 45 km/s. The exponential (left) and Mestel (right) model 
velocities are shown for the cases given in Tables [TJ and [2] i.e. 
Q = 1,1.5,2 as well as A = A m . dx . 

M/L, the curves derived from the Toomre criterion lie at veloci- 
ties that are too high, implying that the picture of a pure disk is 
probably too simple. 

For the case of a very high velocity dispersion, the M/L val- 
ues are substantially higher (columns labelled "cr R = 71 km/s" in 
Tables[TJand[2]i, and now the value from A = A max is also incom- 
patible with the observed stellar range. The velocities show the 
same relative behaviour as for the lower velocity dispersion: the 
model velocities from the Toomre criterion are higher than those 
from A = A lmx , which agree with the measured curve. However, 
the high M/L values might indicate that the actual velocity dis- 
persion is indeed lower than in this extreme case. Of course one 
could argue th at a signific ant amount of dark matter is present as 
well; and vet. iGeha et al.l d2002l) conclude that dark matter does 
not contribute significantly within the half-light radius for their 
(admittedly small) Virgo dE sample. 

The M/L values from A = /l max are thus the only ones that are 
at least partly consistent with the inferred stellar value. A major 
difference between the Toomre criterion and the density wave 
approach is that velocity dispersion does not enter the latter di- 
rectly - only implicitly through the asymmetric drift correction 
- but it does enter the Toomre parameter directly, which yields 
M/L values that are too high. This might hint at a significant dy- 
namically hot component being present in addition to the disk. 
If such a component contributed substantially to the measured 
velocity dispersion, the disk component's own velocity disper- 
sion would be significantly smaller. This would reduce the sur- 
face mass density inferred from the Toomre criterion towards the 
value derived from A = A mdx . 

As for the rotational velocity, since the contribution of the 
asymmetric drift is directly related to the disk's velocity disper- 
sion, its contribution would be much smaller in the case just de- 
scribed. This would lower both the observational and model ve- 
locities, again with particular relevance to the Toomre criterion. 

We illustrate this possibility by repeating our analysis with 
a significantly lower velocity dispersion of cr R = 10 km/s and a 
smaller half-light radius of 0.9 kpc. The latter represents only the 
"inner" component of a double-exponential fit to the light profile 
of Fig. [TJ which yielded a ratio of the two components' scale- 
lengths of 3:10, with a central intensity ratio of 33:2. Obviously, 
the asymmetric drift contribution is much smaller in this sce- 
nario. The choice of cr R = 10 km/s is such that the two ap- 
proaches, A = A mdx and the Toomre criterion, agree with each 



other in the predicted rotational velocities and also the mass-to- 
light ratios (Tables Q]and[2] rightmost column)Q In this case, the 
derived M/L values appear to be somewhat lower than the in- 
ferred stellar M/L — but, given that now the disk contributes 
only part of the light that is used for the calculation, the values 
are very compatible. The model rotational velocities also agree 
with the measured value. This leads to a coherent picture, in 
which IC 3328 consists of a dynamically hot component with an 
embedded thin disk. 

Nevertheless, this interpretation must be regarded with some 
caution, given the measurement uncertainties, the moderate ra- 
dial extent of the kinematical data, and the number of assump- 
tions and simplifications made in the analysis. Clearly, it would 
be desirable to obtain kinematical data of higher quality, in or- 
der to decrease the measurement errors and to trace the curve to 
larger radii. Furthermore, minor axis data would be desirable for 
pinning down the velocity anisotropy. 

Finally, we mention again that no thin disk component is 
seen in the vertical light profile of the apparently edge-on dE(di) 
IC 3435 (see Sect. [TJ. It is of course not known whether this 
galaxy hosts similar spiral structure. Nevertheless, if indeed 
IC 3328 had a thin disk surrounded by a hotter component, then 
IC 3435 would not be an edge-on version of the same kind of 
object. 

5. Summary 

We attempted to gain insight into the actual nature of early-type 
dwarf galaxies with weak spiral structure by analysing the kine- 
matical data of IC 3328, the prototype of this galaxy population. 
To our knowledge, no other object of this kind has kinematical 
data of sufficient quality. Based on mass-to-light ratios derived 
with different methods, we find that the observed velocity dis- 
persion of IC 3328 cannot be fully attributed to a stellar disk, but 
that a distinct dynamically hot component is present. An unam- 
biguous conclusion is, however, not possible, due to the moder- 
ate radial coverage and the lack of minor-axis kinematical data. 
We thus see our study as motivation to increase the amount of 
good-quality kinematical observations for dEs. 
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Appendix A: Finite disk thickness correction of the 
observed rotation curve 

To obtain a first-order correction for the finite thickness of the 
disk, we approximate the density profile of the galaxy with a 
double exponential, 



v = e- R / h * e - z/h * 



(A.l) 



and then integrate along the line of sight over the density- 
weighted rotational velocity: 



1 f + °° 

fi.o.s. = — I dlv cos(ei. . s ,,e^)t' c , 

*J -co 



(A.2) 



with N the normalization factor, R 



— sjx 2 + y 2 , 



and I 



~<Jy 2 + z 2 , such that the x-direction is perpendicular to the line 
of sight. The unit vectors ei. o s . and e,p are defined in the direc- 
tion of the line of sight and the direction of the (cylindrically 
assumed) circular rotation of the galaxy, respectively. 

We assume that the vertical scalelength h z is much shorter 
than the radial scalelength /jr, and perform a Taylor expansion of 
the integrand about z — to second order. After the integration, 
we neglect terms higher than second order in h z , yielding 



v r sin i 



1 r- tan i 

R 1 



d(\n v c ) 
d(lnR) 



(A3) 



Adopting /ir//i z = 3 dLisker et al.ll2006bh and using the inclina- 
tion i = 30° (Sect. 12. Il l, this results in 



fi.o.s. = sin i 



l K 
1-— 4 li 

27 R 2 



d(ln v c ) 
d(lnR) 



(A.4) 
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Fig. B.I. The Toomre Q parameter corrected for the finite thick- 
ness /i z Merit of a galactic disk. The right dashed line repre- 
sents the case <x R = 45km/s, and the left dashed line is for 
<xr = 71 km/s. See the text of Appendix IB1 for details. 



where k denote s the wave number of the Fourier term. Following 
Toomrel d 19641) we consider a disk of finite thickness as a super- 
position of infinitesimally thin disks, 



®(z) 



27rGS k ^ ikx 



f 



dhe 



-k[z-h| ( 



; -|h|/h z 

2h 7 



which leads at the midplane to 
27rG2 k 1 



<D(z = 0) 



1 + kh 7 



(B.2) 



(B.3) 



In the following we use the reductio n term (1 + kh z ) ~ l , which 
is slightly different from that adopted bv lToomrad 19641) . to mod- 
ify the diagnostics of the dynamical state of a galactic disk for 
fini te thickness. 

IFuchs & von LindenI d 19981) show in their Appendix Al how 
such a correction term must be applied to the dispersion re- 
lation for ring-like p erturbations of the disk. In their analysis, 
IFuchs & von LindenI (1 19981 ) adopted an exponential velocity dis- 
tribution. If we swi t ch ba ck to a Schwarzschild distribution as as- 
sumed by iToomrd (1 19641) . we find for the line separating expo- 
nentially growing perturbations from neutrally stable perturba- 
tions in the space spanned by cr-gk^jK (or Q) and the wavenum- 
ber, 



Appendix B: The Toomre stability criterion 
corrected for the finite thickness of galactic 
disks 

The gravitational potential of a Fourier component of the surface 
density of an infinitesimally thin disk is given by 



27rGEk 



ikx-k|z| 



(B.I) 



1 = 



27rGSk 



1 - e 



-<rlk 2 /K : 



1 



1 +kh 7 



(B.4) 



lis 



where Iq denotes the modified Bessel function. Equation [ 
a modified version of equation (62) of Toomr^ ( 1 19641) . 

When Q is drawn as a function of A — 2n/k, we always find 
a lower limit for Q, above which perturbations are stable on all 
wavelengths. This is illustrated as a solid line in Fig. IB. II As 
can be seen from the figure, the allowed Q values are lower than 
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one, which means that a disk of finite thickness can be signifi- 
cantly more massive than an infinitesimally thin disk, but still be 
dynamically stable. However, it needs to be taken into account 
that the vertical scale height h z in units of A alt enters into the 
calculation. Therefore, if one increases the disk surface density, 
Q will drop, but h z /A CI i t drops as well. Thus the most massive, 
but still dynamically stable disk model is reached when Q (orig- 
inally 2=1) has dropped by the same fraction as the originally 
calculated /i z //t C iit- These models are indicated by the intersec- 
tions of the solid line and the dashed lines in Fig. IB. 11 

Appendix C: Density wave theory analysis with a 
correction for finite thickness 

We use the same reduction factor as derived in Appendix iBl to 
correct the prediction of density wave theory for finite thick- 
ness of the dis k. Since we follow "swing-amplification" theory 
(Toomre 1981), we can straightforwardly use the formalism of 
Fuchsl d2001l) . We multiplied the kernel 7C of his equation (68) 



Volterra integral equation. The main result is that the operation 
characteristics of the swing-amplifier mechanism, i.e. at which 
k y and A^ ff amplification peaks, are hardly affected by this correc- 
tion, the effect being of the order of 15%. Only the amplification 
amplitude is much reduced as compared to infinitesimally thin 
disks. This shows again that disks of finite thickness are dynam- 
ically more stable than razor-thin disks. 




